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The van der Waals potential of two atoms in the presence of an arbitrary arrangement of dis- 
persing and absorbing magneto-electric bodies is studied. Starting from a polarizable atom placed 
within a given geometry, its interaction with a second polarizable/magnetizable atom is deduced 
from its Casimir-Polder interaction with a weakly polarizable/magnetizable test body. The general 
expressions for the van der Waals potential hence obtained are illustrated by considering first the 
case of two atoms in free space, with special emphasis on the interaction between (i) two polarizable 
atoms and (ii) a polarizable and a magnetizable atom. Furthermore, the influence of magneto-electric 
bodies on the van der Waals interaction is studied in detail for the example of two atoms placed 
near a perfectly reflecting plate or a magneto-electric half space, respectively. 
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I. INTRODUCTION 

The van der Waals (vdW) interaction of two neutral, 
unpolarized, but polarizable atoms is a well-known con- 
sequence of quantum ground-state fluctuations. For suffi- 
ciently small separations, its physical origin may be seen 
in the electrostatic Coulomb interaction of the atoms' 
fluctuating dipolc moments. The vdW interaction was 
first calculated in this nonretarded limit by London on 
the basis of perturbation theory, who found an attractive 
potential proportional to r^^, where r denotes the inter- 
atomic separation [1]. For larger separations, the vacuum 
fiuctuations of the (transverse) electromagnetic field also 
contribute to the interaction. This was first taken into ac- 
count by Casimir and Polder by means of a normal-mode 
expansion of the electromagnetic field, generalizing the 
London potential to arbitrary distances and showing that 
in particular in the retarded limit the potential varies as 
r-' [2]. 

The theory has since been extended with many re- 
spects, and various factors affecting the vdW interaction 
have been studied. It has been shown that in the case 
of one [3] or both atoms [4, 5] being excited, the vdW 
potential varies as and r"^ in the nonretarded and 
retarded limits, respectively. Thermal photons present for 
finite temperature have been found to lead to a change 
of the retarded vdW potential of two ground-state atoms 
from a r~^- to a r~^-dependence as soon as the inter- 
atomic separation exceeds the wavelength of the dom- 
inant photons [6-9]. The influence of external electric 
flelds on the vdW interaction was addressed, where it 
has been found that the resulting potential varies as 
in the nonretarded limit when the applied fleld is unidi- 
rectional [10]. Generalizations of the vdW interaction to 
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the three- [11-14] and iV-atom case [15, 16] were studied 
flrst in the nonretarded limit and later for arbitrary in- 
teratomic separations, where the potentials were seen to 
depend on the relative positions of the atoms in a rather 
complicated way. 

The two-atom vdW interaction may be strongly af- 
fected by the presence of magneto-electric bodies. This 
was flrst demonstrated by Mahanty and Ninham, who 
employed a semiclassical approach to obtain a general ex- 
pression for the vdW potential of two ground-state atoms 
in the presence of electric bodies [17-19], and applied it 
to the case of two atoms placed between two perfectly 
conducting plates [18]. The situation of two atoms be- 
tween two perfectly conducting plates was later reconsid- 
ered taking into account flnite temperature effects [20]. 
Other scenarios such as two atoms placed within a pla- 
nar electric three-layer geometry [21] or two anisotropic 
molecules in front of an electric half space or within a 
planar electric cavity have also been studied [22]. 

Bearing in mind that the vdW potential of two polar- 
izable atoms may be modified due to finite temperature, 
external fields, or the presence of electric bodies, but re- 
mains attractive in all of these cases, it is rather sur- 
prising that the interaction of a polarizable atom with 
a magnetizable one is repulsive. This was first realized 
by Feinberg and Sucher who restricted their attention to 
the retarded case and found a repulsive potential propor- 
tional to [23]. Their result was later extended to all 
distances [24, 25] , and in particular it was shown that the 
nonretarded potential is proportional to [26]. 

Van der Waals interactions of two atoms exhibiting 
electric as well as magnetic properties have so far only 
been studied in the free-space case. A much richer range 
of phenomena is to be expected when allowing for the 
presence of magneto-electric bodies, where a complex in- 
terplay of the electric and magnetic properties of the 
atoms and the bodies infiuences both sign and func- 
tional dependence of the two-atom vdW potential. This 
problem is addressed in the current work, where we de- 



2 



rive the vdW potential of a polarizable atom with an- 
other polarizable or magnetizable one by starting from 
its Casimir-Polder (CP) interaction with a weakly polar- 
izable or magnetizable body, respectively. This approach, 
which has the advantage of being much simpler than per- 
turbative methods and easily applicable to magnetizable 
atoms, renders general expressions for the two-atom vdW 
potentials of polarizable and/or magnetizable atoms in 
the presence of an arbitrary arrangement of magneto- 
electric bodies, as is shown in Sec. II. In Sec. Ill, the 
general results are applied to the vdW interaction be- 
tween polarizable/magnetizable atoms in free space and 
in front of a magneto-electric plate. A summary is given 
in Sec. IV. 



II. GENERAL THEORY 

Consider first a neutral, nonpolar, ground-state atom 
or molecule A (briefly referred to as atom in the follow- 
ing) in the presence of an arbitrary arrangement of dis- 
persing and absorbing magneto-electric bodies. The atom 
is characterized by its centcr-of-mass position ta and its 
frequency-dependent electric polarizability aA{ijj), while 
the bodies are given by their macroscopic (relative) per- 
mittivity e(r, uj) and permeability /i(r, lu), which are spa- 
tially varying, complex-valued functions of frequency, 
with the corresponding Kramers-Kronig relations being 
satisfied. 

Due to the presence of the bodies, the atom will be sub- 
ject to a CP force F^. Within the framework of macro- 
scopic QED in linear, causal media and by using leading- 
order perturbation theory, it can be shown that this force 
follows from the associated potential [27, 28] 



UA{rA) = ^ r duu^aA{iu)TiG^'\rA,rA,iu) (1) 
according to 



Fa = -VAUA{rA) 



(2) 



(Va = Vr^). In Eq. (1), G^^^r,r',Lu) is the scattering 
part of the classical Green tensor of the electromagnetic 
field, 



G(r,r',Lj) = G(°)(i-,r',t^) + G(i)(r,r', 



(3) 



[G^°\r,r' ,u}), bulk part], which is the solution to the 
equation 



VxK(r,a;)Vx -e{r,aj) 



G{r,r',u;)=S{r-r') (4) 



[K(r, w) = n ^(r, w)] together with the boundary condi- 
tion 



G(r,r',a;)^0 for 



oo. 



(5) 



In order to derive the vdW interaction of atom A 
with a second polarizable atom B, we now introduce 



an additional, weakly polarizable body of (small) vol- 
ume Vb, consisting of a collection of atoms of type B. 
Provided that the atomic number density rjB{v) is suffi- 
ciently small, the electric susceptibility of the additional 
body can be approximated by 



Xe(r,a;) = 



Eq ^r]B{r)aB{oj) if r e Vb, 
iir^Vb, 



(6) 



so that the permittivity of the total arrangement of bod- 
ies reads e(r, w) + Xo(r, uj), and the corresponding Green 
tensor is given by Eq. (4) with e{r,Lu) +Xc(r,Ci;) instead 
of e(r, u)). A linear expansion of this differential equation 
in terms of Xe reveals that the presence of the additional 
body leads to a change of the Green tensor, whose leading 
term is 



AG(r, r',Lu) 

= J d3sXe(s,a;)G(r,s,a;)-G(s,r',a;). (7) 



According to Eq. (1), the resulting change of the CP po- 
tential is 



AUAirA) 



du (^—^ aAiiu) I d'^sXo(s,«u) 

(8) 



27reo Jo 

X Tr[G(rA,s,zu)-G(s, rA.m)]. 



Recalling Eq. (6), on can easily see that AUa is just an 
integral over two-atom potentials C/AB(rA, i"s), 

AUA{rA) = f d^rBVB{rB)UAB{rA,rB), (9) 



where 



UAB{rA,r-B) 



aA{iu)aB{iu) 

X Tl[G{TA,TB,iu)-G{vB,VA,iu)\. (10) 



is the vdW potential between two polarizable atoms in 

the presence of an arbitrary arrangement of dispersing 
and absorbing magneto-electric bodies. The total force 
acting on atom A{B) due to atom B{A) and the bodies 
is just the sum of the single-atom CP force (2) and the 
two-atom vdW force 



Fab(ba) — -'^a(b)Uab{ya,yb), 



(11) 



where in general ¥ab^ —Fba, due to the presence of the 
bodies. Equation (10) agrees with the result obtained on 
the basis of fourth-order perturbation theory, with the 
perturbative calculation being much more lengthy [29, 
30] . Needless to say that the method presented here can 
be easily extended to derive A^-atom potentials [29, 31], 
whereas the perturbative method becomes increasingly 
cumbersome for large N . 

The vdW potential between a polarizable atom A and a 
magnetizable atom B can be derived in a very analogous 
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way. Again wc start from atom A placed within an ar- 
bitrary arrangement of magneto-electric bodies, Eqs. (1) 
and (4), but this time wo add a weakly magnetizable 
body, consisting of a collection of magnetizable atoms of 
type B. For sufBciently small atomic number density the 
magnetic susceptibility of this body approximately reads 

with (3b{^) denoting the magnctizability. The inverse 
permeability of the total arrangment of bodies reads 
K(r,w) — Xni(r,(jj), so that the Green tensor corre- 
sponding to this arrangement is given by Eq. (4) with 
k{t,lo) — x^{y,uo) inwtcad of K{r,w). A linear expansion 
in terms of Xm(r, w) leads to 



AC?(r,r',a;) = 

d^s Xin(s, w) [G(r, s, w) x Vs] • Vs x G(s, r', w). 

(13) 



Combining this with Eq. (1), we find that the change in 
the CP potential due to the presence of the additional 
magnetizable body is 

MJAivA) = J UAiiu) yd^sxm(s,iu) 

X Tr|[G(rA,s,m)x Vs]-VsXG(s,rA,m)|. (14) 

Finally, upon using Eq. (12), Eq. (14) can be rewritten 
in the form of Eq. (9), where now 

UAB{rA,rB) = [ duf-) a^(m)/3B(«M) 
27reo Jo ^ c / 

X Tr[G(r^,r,m)x Vr-VrXG(r,rA,m)]^^^^ (15) 

is the vdW potential between a polarizable atom A and a 
magnetizable atom B in the presence of an arbitrary ar- 
rangement of magneto-electric bodies. To our knowledge, 
Eq. (15) has never been derived so far. 



III. EXAMPLES 

A. Free space 

In order to illustrate the two-atom vdW potentials (10) 
and (15), let us first consider two atoms in free space, 
with the Green tensor being given by [32] 

p—up/c 

G(")(r,r',z«) = -^{a[c/{up)]I-b[c/iup)]epep}, 

''^ (16) 
a{x) = l + x + x'^, b{x) = 1 + 3x + Sx"^ (17) 



(p = r — r'; p = \p\: Cp = p/p; I, unit tensor). In this 
case, the vdW potential between two polarizable atoms, 
Eq. (10), reads 



U^AliT^A,rB) = - 



h_ 

32^ 



2yg / duaA{iu)aB{iu)g{ul/c), 
0' Jo 

(18) 

(x) =2e~^''{S + 6x + 5x^ + 2x^ + x^) (19) 



(Z= |rA — rs|), in agreement with the well-known result 
of Casimir and Polder [2] . From 



g{ur/c) 



4f , 



^[e-^^(9 + 18a;+ 16a;^ 



+ Sa;^ + 3a;* + x^)] 



x=ur/c 



(20) 



it can be seen that the potential (18) is always attractive. 

In the retarded limit, where Z;:g>c/a;min (wmin denoting 
the minimum of all resonance frequencies of atoms A and 
B) the exponential factor effectively limits the u-integral 
in Eq. (18) to a region where 



aA{iu) ~ aA(0), aB{iu) ~ aB(0), 



(21) 



hence the integral can be performed in closed form to 
yield 



U\^{va,Yb) = -jf, 



G7 



2'ihcaA{i))aB{0) 



(22) 
(23) 



In the opposite, nonretarded limit, where I <C c/wmax 
('^max denoting the maximum of all resonance frequencies 
of atoms A and B), the factors aA{iu) and aB{iu) effec- 
tively limit the u-integral in Eq. (18) to a region where 



g{ul/c) ~ g{0) = 6, 
so that the London potential [1] is recovered, 

rr(0) / \ ^6 

U\>j{YA,rB) = —j^ , 



'AB 



Jo 



du aA{iu)aB{iu). 



(24) 

(25) 
(26) 



To calculate the potential between a polarizable atom 
and a magnetizable one, we first use Eq. (16) to derive 



VxG(°)(r,r',m) 



A-Kp^ ("'"''' c ) 
«<».,..'...,XV'.£|(1.^) 



epX/, 



Jxcp. 



(27) 
(28) 



Substituting Eqs. (27) and (28) into Eq. (15) and using 
the identity 



Tr[epX/xep] = —2, 



(29) 
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we find that 



327r3£oZ^ Jo 



aA{iu)PB{iu)h{ul/c), (30) 



/i(x) = 2e-^^(l + 2x + x^), 



(31) 



in agreement witli results found earlier [24, 25]. 

In contrast to tlie attractive vdW potential between 
two polarizable atoms, Eq. (18), tlie vdW potential be- 
tween a polarizable and a magnetizable atom, Eq. (30), 
is always repulsive, as can be seen from 



h{ur/c) 



4f 



— [e-^^(2 + 4a; + 3a;^ +a;^ 



In particular, one finds, on using Eqs. (21) and 

h{ul/c) ~ h{0) = 2, 
that Eq. (30) reduces to 

7t(o) /- ^ ^7 



J x=ur/c' 

(32) 
(33) 



C7 = 



7flc/xoaA(0)/3B(0) 
647r3£o 



and 



t^i°l(rA,rB) 



/4 ' 



167r''eo 



— j aA{iu)l3B{iu) 



(34) 

(35) 

(36) 
(37) 



the nonzero elements of the scattering Green tensor are 
given by (App. A) 



G^%y^rA,rB,^u) = ^J^ dgge"-- 



1 



Jo(gX)(±)J2(gX) b[Jo(gX)(+)J2(gX)] ■ 

b ""^ fc2 



1 

^izix) (r A , rs , m) = (+) — dq q 
G^zz{^A,YB,iu) = -— dqq- 



(38) 



2 -bz+:AM 

^2 



(39) 



[J,y(.T), Bessel function; Z^ = za + zb', X=xb — xa\, where 
ra = ra{q,u) (a = s, p) are the reflection coefficients of the 
half space [cf. Eqs. (46), (63) and (64) below], and 



b = b{q,u) = y-^ + q"^, 

k = k{q, u) = — . 



(41) 
(42) 



According to the decomposition (3) of the Green ten- 
sor, the two-atom potential (10) can be split into three 
parts. 



UABivA, Vb) = UflivA, Vb) + U'J^^iv A, Vb) 
+ Uf^{vA,VB), 



(1), 



(43) 



in the retarded and nonrctardcd limits, respectively. The 
vdW potential between a polarizable atom and a mag- 
netizable one hence shows a power law in the non- 
rctardcd limit which is weaker than the corresponding 
vdW potential of two polarizable atoms by a fac- 
tor of 7/23; while in the nonretarded limit, the potential 
between a polarizable and a magnetizable atom follows 
a power law which is more weakly diverging than 
the corresponding potential between two polarizable 
atoms. 



B. Semi-inflnite half space 

Let us now study the influence of the presence of 
magneto-electric bodies on the two-atom vdW potential. 

To that end, we consider two polarizable atoms placed 
near a homogeneous semi-infinite half space. We choose 
the coordinate system such that the z axis is perpendic- 
ular to the plate with the origin being on its surface, and 
the two atoms lie in the xz plane (Fig. 1). In this case. 



where the bulk-part contribution U^2^{va,vb) is simply 
the free-space result (18), 



J7y2(rA,rB) = --^ / dw (-) aA{iu)aB{iu) 

X Tr [G(0) (r^ , rs , m) • C?(i) (r^ , r ^ , iu)] 
= - oo 3^2, / dM(-) a^(m)aB(w)e-"'/^ 

X / dqqe-^^+ 
Jo 



2/(0-5(0^ 



I- 



br„ 



b fc2 



Ts brp 

b 



bk^ 

J2{qX) 



MqX) 



(44) 



comes from the cross term of bulk and scattering parts 
[with Z = zb — z^, ^ = c/(Zw), and a{x) and b{x) as given 
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X 



A Z 



z = Q 



FIG. 1: Two atoms near a perfectly reflecting plate. 



in Eq. (17)], and 



U^ab{ya,yb) 



— ^ / dw f-) aA{iu)aB{iu) 
SttEq Jo ^ c / 

X Tr [g(i) {va ,VB,iu)-G^^^ (rs , va , iu)] 

fi f°° f°° 



bb' ^ fc4 i + l^lj, I ^^2 



b'k"^ 



MqX)Jo{q'X) + P Ji(gX)Ji(g^X) 



+ 



rjr^ bb'rpr'p b'r.r'p bvpr', 
bb' k'^ bk'^ b'k"^ 



J2{qX)J2iq'X 

(45) 

is the scattcring-part contribution [b' = b{q',u), r'^ = 
ra{q',u) for a = s,p]. 



1. Perfectly reflecting plate 

For a perfectly reflecting plate, the reflection coeffi- 
cients are simply given by 



(46) 



where the upper (lower) sign corresponds to a perfectly 
conducting (permeable) plate. In the retarded limit, 

where 1,za,Zb > c/wmin, C^is is given by Eq. (22), 

whereas u'^^b [Eq- (44)] and U^2b [Eq- (45)] can be given 

in closed form only in some special cases. If X ^ Z+ 
(cf. Fig. 1), we derive, on using the relevant elements 
of the scattering Green tensor as given in App. A 
[Eqs. (AlO) and (All)], 



U 



(1) 

AB 



± — 



32 (X2 + 6;2)C7 



23 l^^Z+il + Z+f 



(47) 
(48) 



where Cy is given by Eq. (23). Thus, recalling Eq. (22), 
the two-atom vdW potential (43) reads 



Cr 32 jX^ + &f)C, 



z+r 



9l 



^- (49) 



In particular, if 2:^/23 ^1, Eqs. (47) and (48) imply that 



(1) 



U 



AB 
(2) 



AB 



U 



(0) 
AB ' 



(50) 
(51) 



so the presence of the perfectly reflecting plate leads to 
an enhancement of the interaction potential, 



Uab = 



AB' 




(52) 



AB 



for a perfectly conducting or permeable plate, respec- 
tively. 

Quite generally, since the bulk part [first term 

on the r.h.s. of Eq. (49)] is negative, the interaction po- 
tential is enhanced (reduced) by the plate if the seat- 
on the 



b'rgr'p tering part u'^g + U^^ [second and third tc 



(2) 



r.h.s. of Eq. (49)] is negative (positive). In the case of a 
perfectly conducting plate, it is seen that especially for 



Z = 0, briefly referred to as the parallel case, U^2b~^^ 



r(2) 

AB 



is positive, and hence the interaction potential is reduced 
by the plate, whereas for X = 0, briefly referred to as the 



vertical case, t^AB + ^As positive and the interaction 
potential is reduced iff 



r(2) 



zb/za < 4.90, 



(53) 



where, without loss of generality, atom A is assumed to 
be closer to the plate than atom B. It is apparent from 
Eq. (49) that for a perfectly permeable plate U^^+U^^ 
is always negative, and hence the interaction potential is 
always enhanced by the plate. 

In the nonretarded limit, where 1,za,Zb <C c/wmax, 
is given by Eq. (25), and from Eqs. (44) and (45) 
we derive, on making use of the relevant elements of the 
scattering Green tensor as given in App. A [Eqs. (A14)- 
(A17)], 



'-'ab ~ 



[4X'^ - 2Z^Zl + X^{Zl + ^2)] Cg 



3/5/5 



u 



(2) 
AB 



'1% 



(54) 
(55) 



{1+ = ^X^ + Zl), where Cg is given by Eq. (26). 

Hence, the interaction potential (43), reads, on recalling 
Eq. (25), 



Ual 



2Z^Zl+X\Zl + Z^)\Ce, 
3/5/5 



9l 

'+ 
(56) 
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conducting plate 


permeable plate 


parallel case 


+ 




vertical case 




+ 



TABLE I: Sign of U^g for a perfectly reflecting plate. 



Let us again consider the effect of the plate on the 
interaction potential for the parallel and vertical cases. 
In the parallel case, Eq. (56) takes the form 



Uai 



which in the on-surface limit ■ 



approaches 



(57) 



Uai 



(58) 



for a perfectly conducting or permeable plate, respec- 
tively. It can easily be seen that the term U^^ [second 

(2) 

term on the r.h.s. of Eq. (57)] dominates the term U^^ 

[third term on the r.h.s. of Eq. (57)], so C^ab + ^1b P^^" 
itive (negative) for a perfectly conducting (permeable) 
plate, and hence the interaction potential is reduced (en- 
hanced) due to the presence of the plate. 

In the vertical case, from Eq. (56) the interaction po- 
tential is obtained to be 



jj _ Ce 2Ce 



Ce_ 

zl 



(59) 



It is obvious that U^g + U)^g [second and third terms 
on the r.h.s. of Eq. (59)] is negative when the plate is 
perfectly conducting, thereby enhancing the interaction 

r(0) 



r(2) 



potential since [first term in Eq. (59)] is negative, 

of a ] 

positive iff 



In the case of a perfectly permeable plate, U^^ + U^^ is 



££<! + 

ZA (3/2)1/3 



1 



~ 14.82, 



(60) 



where atom A is again assumed to be closer to the plate 
than atom B. 

The enhancing/ reducing effect of the perfectly reflect- 
ing plate on the two-atom vdW potential in the various 

cases considered can be systematized in a simple way. 

Since uf^ and [/^^ are negative in all the cases, the en- 
hancement or reduction of the vdW potential due to the 
presence of the plate depends only on the sign of U^g 
and its magnitude compared to that of Uj^g. Moreover, 

the results for the non-retarded limit (the sign of u'^^ 
being summarized in Tab. I) can be explained by using 

the method of image charges, where the two-atom vdW 
interaction is regarded as being due to the interactions 



between fluctuating dipoles A and B and their images A' 

and B' in the plate, with 



-f^int = Vab + Vab' + Vb/ 



(61) 



being the corresponding interaction Hamiltonian. Here, 
Vab denotes the direct interaction between dipole A and 
dipole B, while Vab' and Vba' denote the indirect inter- 
action between each dipole and the image induced by the 
other one in the plate. According to this approach, the 
vdW potential Uab can be identified with the second- 
order energy shift 



AEab 



(OAKOs|-ffi„t|nA)|ms) 



X (nA|(mB|i?mt|OA)|Os). (62) 



(i52(B), \nA(B)}': atomic eigenenergies and eigenstates, re- 
spectively). In this approach, U^^ corresponds to the 

product of two direct interactions, so it is negative in 
agreement with Eq. (56) , because of the minus sign on the 

r.h.s. of Eq. (62). Accordingly, u'^g is due to the prod- 
uct of two indirect interactions and is also negative -in 
agreement with Eq. (56). The terms containing one di- 
rect and one indirect interaction are contained in U^^ 
and determine its sign. We can hence predict the sign of 
from a graphical construction of the image charges, 
as sketched in Figs. 2-5. 

Figure 2 shows two electric dipoles in front of a per- 
fectly conducting plate in the parallel case. The configu- 
ration of dipoles and images indicates repulsion between 
dipole A{B) and dipole B'{A'), so is positive, in 

agreement with Tab. I. On the contrary, in the verti- 
cal case from Fig. 3 attraction is indicated, i.e., negative 
Uab' which is also in agreement with Tab. I. 

The case of two electric dipoles in front of a perfectly 
permeable plate can be treated by considering two mag- 
netic dipoles in front of a perfectly conducting plate, as 
the two situations are equivalent due to the duality be- 
tween electric and magnetic fields in the absence of free 
charges or currents. From Figs. 4 (parallel case) and 5 
(vertical case) it is apparent that the interaction between 
dipole A{B) and dipole B'{A') is attractive in the parallel 
case and repulsive in the vertical case, again confirming 
the sign of as given in Tab. I. 

When the dipole-dipole separation in Fig. 5 is suffi- 
ciently small compared with the dipole-surfacc separa- 
tions, then the direct interaction between the two dipoles 
is expected to be stronger than their indirect interaction 
via the image dipoles. As a result, U^l will be the domi- 



nant term in u'^^+U^^ and U)^^+U^^b becomes positive. 
However, when the dipole-dipole separation exceeds the 
dipole-surface separations, then the indirect interaction 

(2) 

may become comparable to the direct one, and U\g may 
be the dominant term, leading to negative U^ab~^U^2b- 



(1) , rr(2) 
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A' 



FIG. 2: Two electric dipoles near a perfectly conducting plate 
(parallel case). 



FIG. 4: Two magnetic dipoles near a perfectly conducting 
plate (parallel case). 



B' 



A' 



B 




FIG. 3: Two electric dipoles near a perfectly conducting plate 
(vertical case). 



FIG. 5: Two magnetic dipoles near a perfectly conducting 
plate (vertical case). 



The image dipole model hence gives also a qualitative 
explanation of the condition (60). 



2. Magneto-electric half space 

Let us now abandon the assumption of perfect reflec- 
tivity and consider a magneto-electric half space of per- 
mittivity e(w) and permeability /i(w). In this case, the 
reflection coefficients in Eqs. (44) and (45) are given by 



Ts = ra{q,u) = 
Tp = rp{q,u) = 



li{iu)h — b 



M 



fi{iu)h + 6m ' 



s{iu)b + bu 
where b is deflned by Eq. (41), and 



(63) 

(64) 

(65) 



In the retarded limit, 1,za,zb ^ c/wmin (where Wmin 
now denotes the minimum of all resonance frequencies 
of atoms A and B and the magneto-electric medium) 



we may again approximate the atomic polarizabilities by 
their static values, recall Eq. (21), and similarly we may 
set 



e{iu) ~ £(0), ^,{iu) ~ m(0)- 



(66) 



Replacing the integration variable q in Eq. (44) by 
v = bc/u [cf. Eq. (A9)], one can show that the contri- 
bution U^]^Q to the vdW potential takes the form 



he 



327r3/3£2 



aA(0)as(0) j 



Z^A5_ + {Z^ - 2X2) 



A. 



I 



+ 



P 



+ rA5+ + 1 A. 



14+ -I- ^3+ 



+ 2(t>2-l) 



+ (X2-2Z2)('^ + | 



Z^A5+ + {Z^ - 2X^) 



+ I'^As- + lAi- + As- 



A, 



A 



3+ 



(67) 
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where according to Eqs. (63) and (64), the static reflec- 
tion coefficients are given by 



rs = ra{v) 
Tp = rp{v) 



tiiO)v - ^/e{0)^x{0)-^+v^ 



fi{0)v + ^e(0)/x(0) 
s{0)v - V£(0)m(0) + 



and 



/■OO 

Ak± = Axx^ e-^^ [Jo(Ca;) ± JiiQx)] , 
Jo 

Jo 



(68) 
(69) 

(70) 
(71) 



with X = l + vZ+ and C = Xy/v^ — 1 (for expUcit expres- 
sions of Ah± and Bk, see App. B). Similarly, Eq. (45) 
reduces to 

he 



u 



(2) 



AB 



647r%g 



aA{0)aB{0) I dv J dv' 



X I (rpr'p [3v\'^ - 2{v'^ + v'^) + 2] + r.r', - r.r'pv'^ 

+ {tst's + rpr'pv'^v'^ + r^r'^v'"^ + rpr'^v^)M.^ 
[r^ = ra{v') for (T = s,p], where 

dxa;6e-(''+''')^+"J,(Ca;)J,(C'a;) 



(72) 



(73) 



(C = X^/v'"^ — 1), which can be evaluated analytically 
only in some special cases. In particular, when X <^ 
then approximately 

M,-J,[i))j^ axxe {v + v'yzl' 

(74) 

Analytic expressions for U^g and u'^^ in the nonre- 

tardcd limit, 1.za,zb <C c/[-y/e(0)/i(0) Wniax] [with cjmax 
being the maximum of all resonance frequencies of atoms 
A and B and the magneto-electric medium], can be ob- 
tained by using in Eqs. (44) and (45), respectively, the 
relevant elements of the scattering part of Green tensor 
as given in App. A. In the case of a purely electric half 
space {n= 1) we derive [Eqs. (A20)-(A23)] 

Ce [4X4 - 2^2 ^2 + X^iZ"^ + Zl)]D 
Uab=-^ + 



E 

where Cq is given by Eq. (26), and 



D 



Wn^el JO 



du aA{iu)aBiiu) 



du aA{iu)aB{iu) 



e(iu) — 1 
e{iu) + 1 ' 
e{iu) — 1 



1 2 



e{iu) + 1 



(75) 

(76) 
(77) 



In particular, in the limiting case when 1<^Z+, Eq. (75) 
reduces to 



_ Ce , {X^-2Z^)D 

UAB-—j^+ 



(78) 



It is seen that the second term on the r.h.s. of this equa- 
tion is positive (negative) in the parallel (vertical) case, 
so the vdW potential is reduced (enhanced) by the pres- 
ence of the electric half space. 

In the case of a purely magnetic half space {s=l) we 
derive [Eqs. (A24)-(A27)] 

rr _ Ce . [Z'-2X^ + 3Z+{l+-Z+)]F 
Uab--^ + , (79) 



where 



F = 



647r%2 



(iX \ 2 
— j aA{iu)aB{iu) 



^ [m(w)-1][^(w)-3] 

At(m)-M ■ ^ ' 

(2) 

Note that does not contribute to the asymptotic 

nonretarded vdW potential Uab for the purely magnetic 
half space. In particular in the limiting case when X <^ 
Z+, Eq. (79) reduces to 



Uab = - 



{2Z'- - X'-)F 
2l^Z+ 



(81) 



It is seen that the second term in the r.h.s. of this equa- 
tion is negative (positive) in the parallel (vertical) case, 
so the vdW potential is enhanced (reduced) due to the 
presence of the magnetic half space. 

It should be pointed out that the nonretarded limit 
for the magneto-electric half space is in general incom- 
patible with the limit of perfect reflectivity [e{iu) oo 
or /i(m)^oo] considered in Sec. IIIBl, as is clearly 
seen from the condition given above Eq. (75) [cf. also 
the expansions (A18) and (A19), which are not well- 
behaved in the limit of perfect reflectivity]. As a con- 
sequence, Eq. (79) docs not reduce to Eq. (56) via the 
limit ij{iu) ^oo. It is therefore remarkable that the re- 
sult for a purely electric half space, Eq. (75), does reduce 
to Eq. (56) in the limit e{iu) ^ oo, as already noted in 
Ref. [33] in the case of the single-atom potential. 

Figures 6-8 show the results of an exact (numerical) 
calculation of the vdW interaction between two identi- 
cal atoms near a semi-infinite half space, as given by 
Eqs. (43) together with Eqs. (18), (44), and (45) as well 
as Eqs. (63) and (64). In the calculations, we have used 
single-resonance models for both the polarizability of the 
atoms. 



aA{to) = asito) = lim 



2 ^io|dio|2 



(82) 



(with LOiQ and dio denoting the frequency and electric 
dipole matrix element of the dominant atomic transition. 
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FIG. 6: The normalized vdW potential of two atoms in 

the parallel case in presence of (a) a purely dielectric half 
space with uj-pe/t^io ~ 3, LUTe/^^io = 1, and 7e/a;i() — 0.001 (b) 
a purely magnetic half space with cijpm/<^io = 3, uoTm/L^us = 1, 
and 7m/wio = 0.001 is shown as a function of the atom-atom 
separation I, with the atoms placed at at distance za = zb 
= O.Olc/aiio (solid line), 0.2c/a;io (dashed line), and c/wio 
(dotted line) from the half space. 



FIG. 7: The normalized vdW potential of two atoms in the 

vertical case in the presence of (a) a purely dielectric half 
space and (b) a purely magnetic half space is shown as a 
function of the atom-atom separation /. The distance be- 
tween atom A (which is closer to the surface of the half space 
than atom B) and the surface is za = O.Olc/wio (solid line), 
0.2c/a;io (dashed line), and c/lou) (dotted line). All other pa- 
rameters are the same as in Fig. 6. 



respectively), and the permittivity and permeability of 
the half space, 

£(a;) = 1 + 4^- , (83) 

^i{oJ) = 1 + - — . (84) 

In the figures the potentials and the forces are normalized 
w.r.t. their values in free space as given by Eq. (18), so 
one can clearly see that the vdW interaction is unaffected 

by the presence of the half space for atom half-space sep- 
arations that are much greater than the interatomic sep- 
arations (the curves approaching unity for za,Zb^ l), 
while an asymptotic enhancement or reduction of the in- 
teraction is observed in the opposite limit. 



Figure 6(a) shows the dependence of the normalized 
vdW potential on the atom-atom separation I in the 
parallel case (Z = 0) for different values of the distance 
za = zb of the atoms from a purely dielectric half space. 
The ratio of the interatomic force along the connecting 
line of the two atoms, Fabx{1) [Eq- (H)]) to the corre- 
sponding force in free space, F^bx^^)i follows closely the 
ratio Uab{1)/U^^{1), so that, within the resolution of 
the figures, the curves for Fabk (OZ-^abx (0 (^lot shown) 
would almost coincide with those for Uab{1)IU^ab{^)- ^he 
figure reveals that due to the presence of the dielectric 
half space the attractive vdW potential and force arc re- 
duced, in agreement with the predictions from the non- 
retarded limit, Eq. (78). The relative reductions of the 
potential and the force arc not monotonic, there is a 
value of the atom-atom separation where the reduction 
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FIG. 8: The normalized vdW force acting on atom A (which 

is closer to the surface of the half space than atom B) in the 
presence of (a) a purely dielectric half space and (b) a purely 
magnetic half space is shown as a function of the atom-atom 
separation I. All parameters are the same as in Fig. 7. 

is strongest. 

The Z-dependence of Uab{1)/U^^{1) in the presence 
of a purely magnetic half space in the parallel case 
is shown in Figs. 6(b). The corresponding force ra- 
tio FABxil) / Ffg^{l) (not shown) again behaves hke 

UAB{l)/U^g{l). The figure indicates that the presence 
of a purely magnetic half space enhances the vdW in- 
teraction between the two atoms, with the enhancement 
increasing with the atom-atom separation, in agreement 
with the nonretarded limit, Eq. (81). 

Figure 7 shows C/^b(Z)/C/^°2(0 the vertical case 
{X = 0) when the half space is purely dielectric [Fig. 7(a)] 
or purely magnetic [Fig. 7(b)]. In the figures, atom A is 
assumed to be closer to the surface of the half space than 
atom i?, and the graphs show the variation of the vdW 
potential with the atom-atom separation / for different 
distances za of atom A from the half space. It is seen that 
for a purely dielectric half space the potential is enhanced 
compared to the one observed in the free-space case — in 



agreement with Eq. (78). Note that there are values of 
the atom-atom separation at which the enhancement is 
strongest. 

For a purely magnetic half space, the potential is 
seen to be typically enhanced although for very small 
atom-atom separations a reduction appears [inset in 
Fig. 7(b)] — in agreement with Eq. (81). Due to this slight 
reduction for small atom-atom separations, the relative 
enhancement is not monotonous, in contrast to what is 
suggested by the large figure. 

Whereas the force -Fsyiz (OZ-P'baz (0 for the force act- 
ing on atom B (not shown) again follows closely the po- 
tential ratio Uab{1)/U^ab^^) f*-"" both dielectric and mag- 
netic half spaces (as in Fig. 7, not shown), the ratio 

FABz{l)/F^Bzi^)^ t^i® force acting on atom A notice- 
ably differs from UAB{l)/ufl{l) (Fig. 8). Clearly, the 
difference is due to the fact that the atom A{B) which is 
responsible for the force Fba(ab)z{1) is situated on the 
same side of atom B{A) as the half space in the former 
case, but on a different side in the latter case (cf. Figs. 3 
and 5). 



IV. SUMMARY AND CONCLUSIONS 

Starting from the CP potential of a single polarizable 
atom placed within a given arrangement of magneto- 
electric bodies, we have presented a macroscopic deriva- 
tion of two-atom vdW potentials: By introducing an ad- 
ditional weakly polarizable body and linearly expanding 
the resulting CP interaction in terms of the body's sus- 
ceptibility, the vdW potential between two polarizable 
atoms in the presence of an arbitrary arrangement of dis- 
persing and absorbing magneto-electric bodies has been 
inferred. The vdW potential between a polarizable atom 
and a magnetizable one has been derived in a similar way 
by introducing a weakly magnetizable body. The general 
formulas have been used to study the influence of polar- 
izability and magnetizability on the vdW potential be- 
tween two atoms in free space. In particular, it has been 
shown that the vdW interaction of a polarizable atom 
with a magnetizable one is always repulsive, in contrast 
to the well-known attractive potential between two po- 
larizable atoms. 

To illustrate the influence of the presence of magne- 
todielectric bodies on the vdW potential, we have con- 
sidered the example of two polarizable atoms near a per- 
fectly reflecting plate. It has turned out that due to the 
presence of the plate the attractive vdW interaction be- 
tween the atoms can be enhanced or reduced depend- 
ing on the magneto-electric properties of the plate and 
the specific alignment of the atoms with respect to the 
plate. In particular, in the nonretarded limit these effects 
can be qualitatively explained using the method of im- 
age dipoles. To be more realistic, we have also calculated 
the vdW potential for the case of the two atoms near a 
magneto-electric half space of finite permittivity and per- 
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meability. The analytical results show that in the nonre- 
tarded limit the potential in the case of a purely electric 
half space is reduced (enhanced) compared to its value 
in free space in the case of parallel (vertical) alignment 
of the two atoms, while in the case of a purely magnetic 
half space it is enhanced (reduced) for parallel (vertical) 
alignment of the two atoms. The numerical computation 
of the potential in the whole distance regime confirms the 
analytical results. In addition, it shows that the relative 
enhancement / reduction of the vdW interaction is not al- 
ways monotonous, but may in general display maxima or 
minima, in particular in the case of a purely dielectric 
half space. 

In conclusion, the examples studied in this work sug- 
gest that the sign of the vdW potential is entirely deter- 
mined by the electric/magnetic nature of the interacting 
atoms, while the strength of the respective attractive or 
repulsive potentials can be controlled by the presence of 
magneto-electric bodies. 
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APPENDIX A: SCATTERING GREEN TENSOR 
FOR A SEMI-INFINITE HALF SPACE 

The scattering Green tensor for a semi-infinite half 
space can be given in the form [34] 

G^^\r,r',iu) = J d^ge^i-^'— ^')G(i)(q,z,z',m) (Al) 
(q _L Bz), where 

G(i) (q, z,z',iu) = -^^Yl e-r,e-''(^+^') , (A2) 
with 



a=s,p 



e, = sin (t>ex — cos 6 



(A3) 



tq 



ep =T^{cos (f>ex+ sin (f>ey) - ^e;, (A4) 



{eq = cosfpex + siiKpey = q/q, q ~ \q\) denoting the po- 
larization vectors for s- and p-polarized waves propa- 
gating in the positive(-|-)/negative(— ) z-direction. Fur- 
ther, b and k are defined according to Eqs. (41) and 
(42) , respectively, and the reflection coefficients Vg , Vp 
are given by Eq. (46) for a perfectly reflecting plate and 



by Eqs. (63) and (64) for a magneto-electric half space. 
Equations (A3) and (A4) imply that 



^ sm 

ej'e^ = I — sm0cos0 cos 




(A5) 



-|t COS^ (/> 

I 2 

— -p- sin ^ cos ( 



— p- sm (f> cos ( 

l2 o 

-frsm (p 



ibq 



■ sm ( 



(A6) 



Substituting these results into Eqs. (Al) and (A2), per- 
forming the (/(-integrals by means of [35] 



27r 



'cosivx) = 2Tri''J^{x), (A7) 



and using the relation 

Jl{x) Jq{x) - J2{x) 



(A8) 



we arrive at the Eqs. (38)-(40). 

In the particular case of a perfectly reflecting plate in 

the retarded limit, it is convenient to replace the integra- 
tion variable q in Eqs. (38) (40) in favor oiv = bc/u, i.e., 
q = \/tfl — \u/c [see Eq. (41)], and hence 



f 

Jo 



Aq- 



i: 



dv 



(A9) 



For X <C , the exponential terms effectively limits the 
integrals in Eqs. (38) (40) to the region where qX <^1, 
hence we can approximate J^^qX) by Ji/(0) = 5i/O5 such 
that the nonzero scattering- Green tensor components 
read 



<^l^i (r A , r s , iu) = G^,^^ (r ^ , r s , m) 
1 

1 



87rZ+ 



uZ+ 



I f c 



+ ■ 



2'kZ+\uZ+ v?Z% 



-uZ^ /c 



-uZ-\- / c 

(AlO) 
(All) 



leading to Eqs. (47) and (48), recall Eq. (21). 

In the nonretarded limit it can be shown that the main 
contribution to the frequency integrals comes from the 
region where u/{cb) <C 1 (cf. Ref. [36]). In this region we 
have 



6a 1 



b. 



(A12) 
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By changing the integration variable q according to 
9 



dg- 



roo 

J u/c 



db ... 



(A13) 



and setting the lower limit of integration to zero, from 
Eqs. (38) (40) we find, after some algebra, the nonzero 
elements of the scattering Green tensor to be approxi- 
mately given by 



2X2 



Zl c^r^ 



_ ZXZ+ c^Tp 
2Zl c^r^ 



(A14) 
(A15) 
(A16) 
(A17) 



with h 



X2 + Zl, leading to Eqs. (54) and (55). 



'-+ — Y 

For a magneto-electric half space in the nonretarded 
limit, we apply a similar procedure as below Eq. (All) 
and expand the reflection coefficients given by Eqs. (63) 

and (64) in terms of u/(bc), 



fi{iu) — 1 ij.{iu)[e{iu)fj,{iu) — 1] 
fi{iu) + 1 
e{iu) — 1 



[/i(n/) + 1]2 

e{ru)[:{ia)i.i(iu) - 1] 



s{iu) + 1 



[s{iu) + lY 



(A18) 
(A19) 



Substituting (A18) and (A19) into Eqs. (38) (40) and 
keeping only the leading-order terms of u/bc, in the case 
of the purely electric half space we can ignore rg and the 

second term in the r.h.s. of Eq. (A19), so the relevant el- 
ements of the scattering Green tensor arc approximately 



<^ix(rA,rB,m) = 



2X^ - Z\ e{iu) - 1 
477/^ V? e{iu) + 1 ' 
1 e{iu) - 1 



G 



(1) 



(r A, rB,m) = 
G^^){rA,rB,iu 



47rZ^ V? s{iu) + 1 ' 
_ 3XZ+ (?_ e{iu) - 1 
47r/^ ^ e{iu) + 1 
X2 - 2Zl c2 e{iu) - 1 



AttI\ V? e{iu) + 1 ■ 



(A20) 
(A21) 
(A22) 
(A23) 



For a purely magnetic half space, the first term on the 
r.h.s. of Eq. (A19) vanishes, so the leading order of u/hc 
is due to the second term as well as the first term on 
the r.h.s. of Eq. (A18), so the nonzero elements of the 
scattering Green tensor can be approximated by 



G^^^{rA,YB,iu) = 



1+ — Z+ ^{iu) — 1 
47rX2 ^i{iu) + l 
Z+U - Z^ 



G^yJ{rA,rB,iu} 



1+-Z+ 



[li{iu) - 1] 



+ 



167rX2 
Z+1+ — Z\ fj,(iu) — 1 



47rX2/+ /i(m) + 1 ' 



(A25) 



G 



(1) 

xz{zx) 



{ta, tb, iu) = (±)^t_^[^(i„) _ 1], (A26) 



G'^J'J{rA,rB,iu) 



167r/+ 



[li{iu) - 1]. 



(A27) 



APPENDIX B: EXPLICIT FORMS OF Ak± AND 
Bk IN EQS. (70) AND (71) 

The integrals in Eqs. (70) and (71) can be performed 
to obtain the following explicit expressions: 



^3+ = 



6A 



A4+ = 



(A2 + C2)' 

6(A3 - 4AC2) 

(A2+C2)5 

6(4A2 - C') 
(A2 + C 



2\ 2 



^4- = 



_ 6(4A4-27A2c2-F4C^) 
(A2 + C2)^ 



^5+ 



30(4A3 -3AC2) 

(A2+C2)^ 



A.- = 



B. = 



30(4A5 - 41A3C2 + isAC") 
(A2 + C2)^ 

3A(2A^ -3^^^) 

(A2 + C2)^ 

_ 3(8A4- 24A2C2 +3C^) 
(A2 + C2)^ 

15A(8A'^ - 40A2C2 + IbC^) 
(A2 + C2)^ 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 
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